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The convergence properties of two perturbative schemes to sum the so-called folded diagrams 
are critically reviewed in this work, with an emphasis on the intruder state problem. The methods 
we study are the approaches of Kuo and co-workers and Lee and Suzuki. The suitability of the two 
schemes for shell-model calculations is discussed. 



I. INTRODUCTION 



Recent works on the theory of the effective interaction for finite nuclei, have focussed on properties like hermiticity 
and the order-by-order convergence of the perturbative expansion for the effective interaction [1-4] . Especially, tech- 
niques to sum up subsets of diagrams to infinite order have received increased interest. Two such schemes are the 
summation of folded diagrams by Kuo and co-workers, exposed in ref. |^] (hereafter referred to as the FD method), 
and the approach proposed by Lee and Suzuki || (hereafter referred to as the LS method). The starting point for 
these iterative schemes is to define the so-called Q-box, which in a diagrammatic language means the sum of all 
valence-linked and irreducible diagrams to a given order in the interaction. The Q-box depends on the so-called 
starting energy, defined as the unperturbed energy of the interacting nucleons. 

The purpose of these methods is to calculate an effective interaction H c s in order to reduce the full shell model 
eigenvalue problem, to one which is tractable within a physically selected space, the so-called model space. One then 
obtains a secular equation 



acting solely within the model space, defined by the projection operator P, which projects out the model-space 
components of the true eigenfunction \*ff). The quantity Ep is the model-space eigenvalue, and the d eigenvalues 
Ep (d being the dimension of the model space) are supposed to reproduce the corresponding exact eigenvalues [^). 
The degrees of freedom not represented by the model space are supposed to be taken care of by the perturbative 
expansion. The excluded degrees of freedom are obtained through use of the projection operator Q, with P + Q = 1 
and PQ = 0. The corresponding eigenvalues are Eq. 

The above techniques have been applied extensively to light nuclei, in particular to the nuclei 18 O and 42 Ca 
where the model space is rather simple, two particles in the sd- or p/-shells, respectively. Here, however, a special 
problem occurs, in the literature called the intruder-state problem. In short, this means that the P- and Q-states 
do not separate. Experimentally, it is rather well established that the first excited + state in both 18 and 42 Ca 
is dominated by Q-state components, i.e. shell model configurations outside the sd- or p/-shell. Then a comparison 
between theory and experiment becomes difficult. Using the FD method to calculate H r, eq. (1) should reproduce 
the experimental levels dominated by the P-state configurations. For the two examples mentioned above this means 
that the first excited + state in both 18 O and 42 Ca should not be accounted for by the theory. On the other hand 
the LS method, for a properly chosen starting energy, should be able to reproduce the d lowest states of e.g., 18 
and 42 Ca, irrespectively of whether or not they are dominated by Q-space components. In view of this, one might 
expect the FD effective interaction to be more appropriate for shell-model calculations than the LS one. An important 
application of such an effective interaction is spectroscopic shell-model analyses of systems with more than two valence 
particles, and then only P-space degrees of freedom are considered. 



PPeffP|*) =E P P\V), 



(1) 
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In several papers, e.g., ref. |2), we have investigated the above mentioned nuclei using both the FD and the LS 
methods. Our results show little difference in the final energy spectra. The question then arises how one should relate 
the energy levels calculated by respectively the FD and LS methods to the experimental levels. 

In the present paper we investigate this problem in more detail and explore the convergence properties of the FD 
and LS methods in the presence of so-called intruder states. In the next section we employ a simple 2x2 model to 
study the difference between calculations with a complete and a limited perturbation expansion. Then in section 3 
we apply the FD and LS methods to the nucleus 18 O, in order to test if differences seen in the exactly solvable model 
pertain to the realistic case as well. Of importance here is the fact that the simple 2x2 model allows us to define 
an exact Q-box, whereas in the realistic cases the Q-box must be approximated to a given order in the interaction. 
Finally, our conclusions are drawn in section 4. 



II. INTRUDER STATES IN AN EXACTLY SOLVABLE MODEL 



It is a well-known fact that the presence of so-called intruder states |5||7j] may lead to divergence of the order-by-order 
pertubative expansion for fleff- 

One way to handle intruder states is that of introducing an enlarged model space which includes intruder config- 
urations. Such calculations do however become prohibitively time-consuming for nuclei heavier than sd-shell nuclei. 
There are also other approaches which aim at overcoming the divergence of the order- by-order expansion of H c g. One 
may, e.g., regroup the perturbative expansion and perform infinite partial summations to obtain convergence. Two 
such possibilities are represented by the summation of the folded diagrams by the FD and LS schemes discussed above. 
To better understand the above convergence arguments and the structure of the folded-diagram and the Lee-Suzuki 
methods we will first demonstrate certain properties of these methods by considering a simple 2x2 model, employed 
e.g. by the authors of ref. in their discussion of perturbative many-body approaches. Moreover, the model can be 
used to demonstrate the connection between intruder states and the convergence of the perturbative expansion. Last 
but not least, this model allows us to define an exact Q-box, which in turn can be used to make statements about 
the convergence of the FD and LS iterative schemes. 

However, first we need to repeat the equations pertinent to the FD and LS methods. The effective interaction for 
the FD method is given by 

H eS = H +]im (2) 



n — >oo 



with 

V 



O) _ \ - 
ff ~ I , 



Q m {v^- 1} } m . 0) 



Here we have defined Q m = The energy uj is the so-called starting energy, defined as the unperturbed energy 

of the interacting nucleons. 

The Lee and Suzuki (LS) expansion for the effective interaction is given as Q 

H cS = H Q + lim R n , (4) 

n — >oc 

with 



Rn = 



ri-1 n-1 n 1 



1 — Ql — ^ Qm JJ Rk 

in— 2 k—n—7n+l 



Q. (5) 



Eqs. (^) and (|J) differ not only in structure, but also in convergence properties. The convergence criterion for 
the method of Lee and Suzuki is related to the choice of the so-called starting energy ui only. The LS expansion 
converges to those eigenvalues which are closest to u>, which means that we should be able even to reproduce the 
Q-space eigenvalues given an appropriate starting energy. This means that if the LS expansion converges for the 0^ 
state of, e.g. 18 0, the structure of the exact wave function does not correspond to that of the two-particle space we 
have chosen. In contrast, the FD-method converges to the states dominated by the model-space components. 

In order to study the importance of intruder states, we let our hamiltonian depend linearly on a strength parameter 
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H = H Q + zH u 

with < z < 1, where the limits z = and z = 1 represent the non- interacting (unperturbed) and fully interacting 
system, respectively. The model is an eigenvalue problem with only two available states, which we label P and Q. 
Below we will let state P represent the model-space eigenvalue whereas state Q represents the eigenvalue of the 
excluded space. The unperturbed solutions to this problem are 

H $ P = ep$p (6) 

and 

#o$Q = £q$q, (7) 

with e P < €q. We label the off-diagonal matrix elements X, while X P = ($p| H x |3>p) and Xq = ($q| Hi |$q). The 
exact eigenvalues problem 



6p + zXp zX 
zX e Q + zXq 



(8) 



yields 

£(z) = i {e P + e Q + zip + zX Q ± (e Q - e P + zX Q - zX P ) (9) 



^ _| 4z 2 X 2 

(eQ-ep+zXQ-zXp) 



The authors of ref. B demonstrated how Brillouin-Wigner and Rayleigh-Schrodingcr (RS) perturbation theories relate 
within the framework of this simple model. An RS expansion for the lowest eigenstate (defining states P and Q as 
the model and excluded spaces, respectively) can be obtained by expanding the lowest eigenvalue as 

Z 2 X 2 Z 3 X 2 (X Q -Xp) z 4 X 2 (X Q -Xp) 2 z 4 X 4 , . 

E = e P + zX P + + - P > + 1 Q P > - + ..., (10) 

ep-e<3 (ep-e Q ) 2 (ep - £q) a (ep - eq) d 

which can be viewed as an effective interaction for state P in which state Q is taken into account to successive orders 
of the perturbation. In this work we choose the parameters ep = 0, eq — 4, Xp = —Xq = 3 and X = 0.2. The 
exact solutions given by eq. (||) are shown in fig. 1 as functions of the strength parameter z. Pertinent to our choice 
of parameters, is that at z > 2/3, the lowest eigenstate is dominated by $q while the upper is <i>p. At z = 1 the <3?p 
mixing of the lowest eigenvalue is 1% while for z < 2/3 we have a $p component of more than 90%. The character of 
the eigenvectors has therefore been interchanged when passing z — 2/3. The value of the parameter X represents the 
strength of the coupling between the model space and the excluded space. Thus, this simple model allows us to study 
how the perturbation expansion with a model space defined to consist of state P only, behaves as the interaction 
strength z increases. The order-by-order convergence in eq. ([!(]) was discussed by the authors of ref. Here we 
will thence only repeat their conclusion: For small values of z one obtains good convergence to the lower eigenvalue. 
For larger values of z (i.e. z > 2/3), increasing orders in the perturbation expansion yield a divergent perturbation 
series, as expected 0. However, it may be possible to rewrite the perturbative expansion in such a way that one sums 
subsets of diagrams to all orders. The hope is then that the expansion becomes convergent for appropriate infinite 
partial summations. This is actually the philosophy behind both the FD method and the LS method. Having defined 
a set of linked and irreducible valence diagrams, the so-called Q-boxQ, we can define an iterative scheme to sum the 
contributions from folded diagrams. The Q-box serves therefore as the starting point for our iterative schemes, and 
within the framework of our 2x2 model wc can study the FD and LS methods for various approximations of the 
Q-box. As an example, to fifth order in the parameter z we have 



A z 2 X 2 z 3 X 2 X Q z 4 X 2 X 2 z 5 X 2 X Q 

Q b = zX P + + - % + ? % + 7 ^, (11) 

ep - £Q (ep - tQ) (ep - hq) 6 (ep - £q) 4 



and it is easy to see that a Q-box of order I + 2 can be written as 

z 2 X 2 




which in the limit I — > oo givesQ 



The Q-box should not be confused with the exclusion operator Q. 
2 1 — gives a second-order Q-box, I = 1 a third-order Q-box and so forth. 
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z 2 X 2 

<9cxact(ep) = zX P + — , (12) 

ep - £q - zXq 



if 

zXr 



< 1. (13) 

ep - £Q 

The latter equation clearly restricts the possible values of ep for given Xq and eg. Actually, if we let ep vary, our 
choices for Xq and eg restrict ep to ep < 1 and ep > 7, in order to have a finite Q-box. 

The exact Q-box may be used to define the first iteration of the FD expansion as a function of the starting energy 
to as^| 

Ai = UJ + QcxactM (14) 



The subsequent steps are 



and 



In general we have 



A 2 = Ai + ^ ~-^{\ x -u) m , (15) 



m! duj r 



. _ 1 A rn Cl 

+ (16) 

m— 1 
°° 1 r} m (l 

= + £^(^-») m - (17) 

m— 1 

If this iteration scheme converges we have, with A = A„ = A„_i, 

A = cj + zXp + - — , (18) 

A - eg - zX Q 

which is just eq. ([)]), so that we have obtained the true eigenvalues. In a similar way we can use the LS expansion 
defined in eq. ([!]) to sum the folded diagrams. 

TABLE I. The exact solutions Ep ( model space) and Eq ( excluded space) of eq. (9) as functions of the strength parameter 
z. The results obtained with the LS and FD methods with an exact Q-box as functions of the starting energy u) are also given. 





z 


E P 


Eq 


FD 


LS 


0.5 


0.0 


0.00 


4.00 


0.00 


0.00 




0.2 


0.60 


3.40 


0.60 


0.60 




0.4 


1.20 


2.80 


1.20 


1.20 




0.6 


1.77 


2.23 


1.77 


1.77 




0.7 


2.17 


1.83 


2.17 


1.83 




0.8 


2.43 


1.57 


2.43 


1.57 




1.0 


3.02 


0.98 


3.02 


0.98 


7.5 


0.0 


0.00 


4.00 


0.00 


4.00 




0.2 


0.60 


3.40 


0.60 


3.40 




0.4 


1.20 


2.80 


1.20 


2.80 




0.6 


1.77 


2.23 


1.77 


2.23 




0.7 


2.17 


1.83 


2.17 


2.17 




0.8 


2.43 


1.57 


2.43 


2.43 




1.0 


3.02 


0.98 


3.02 


3.02 





3 We replace here ep with uj. 
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We demonstrate the properties of the LS and FD methods with an exact Q-box in table 1 for two starting energies, 
0.5 and 7.5 (arbitrary units). Clearly, we see that with a starting energy 0.5 the LS method yields the Q-space 
eigenvalue at z > 2/3. Below z — 2/3, the LS method reproduces the P-space eigenvalue. With a starting energy 
of 7.5, the LS method gives the Q-space eigenvalue for z < 2/3, whereas the P-space eigenvalue is reproduced for 
z > 2/3. Thus, this simple model demonstrates nicely the properties of the LS scheme, i.e. it converges to those 
eigenvalues which are closest to the chosen starting energy, irrespectively of the structure of the wave function. From 
table 1 we also see that the FD method always reproduces the P-space eigenvalue. At z = 2/3, the FD scheme does 
not stabilize, the eigenvalue fluctuates and there is no convergence. More precisely this means that we can not go from 
the lower to the upper eigenvalue along the real axis z [|| . The fluctuation is intimately connected to the convergence 
criterion for the FD scheme. If one of the solutions contains more than 50% valence state intensity, the FD method 
converges to that solution. At z = 2/3, equal admixtures of $p and <&q occur in the true wave functions. 

This simple example, starting with the exact Q-box, serves to demonstrate significant differences in convergence 
behavior of the two methods. Of importance is the fact that with the LS scheme we are able to reproduce a Q-space 
state insofar we define a starting energy which is close to the actual Q-space state. However, in actual nuclear structure 
calculations, we are not able to define an exact Q-box. The question we wish to bring to the attention in this work, 
is if an approximation of the Q-box still gives the same difference between the LS and FD methods at z = 1. To 
shed light on this, we exhibit in fig. 1 results obtained with various low-order approximations to the Q-box, recall eq. 
([TT]) . The starting energy for the LS and FD calculations was set equal 0.5, since we are interested in seeing if the 
convergence criterion for the LS method holds for an approximate Q-box. Clearly, with a Q-box of second order in 
the interaction, the difference^ between the LS and the FD method is negligible. A second-order Q-box corresponds 
to setting Xq = 0, which means that in the range of physically interesting z values (0 < z < 1) there is no intruder 
state, yielding an approximately straight line for the LS method. Up to fourth order in the interaction, both the LS 
and FD methods yield almost the same value. With a fifth-order Q-box the FD method becomes unstable at z = 1. 
This can be inferred from the structure of the FD expansion, which with a fifth-order Q-box reads 

z 2 X 2 z 3 X 2 X Q ^X 2 X 2 z 5 X 2 X* 
A„ = e P + zX P + - + + ->- 



K-l - £Q (K-l - £q) 2 (A„-i - eg) 3 (A„_i - £q) 4 

The first iteration is just Ai = ep + Q, which gives a result close to the model space eigenvalue « 3 at z = 1. Inserting 
Ai and higher iterations in the above expansion results in an FD expansion which for each iteration will fluctuate 
between a large value A and ep + Xp. The large eigenvalue stems from the fact that the energy denominators in 
the higher-order terms become small. With X — 0.2 this oscillating behavior sets inn already at fifth order in the 
interaction, whereas if we choose the coupling between the model space and the excluded space to be X — 0.01, 
this divergence appears first with a Q-box of tenth order. However, we will ultimately end up with a series which 
fluctuates, except for the trivial case X — 0. With the given parameters and a Q-box of tenth or higher order in 
the interaction, the FD method converges only if z < 2/3. If the FD method converges, we obtain a result close to 
the model-space eigenvalue, and, if the coupling between the model space and the excluded space is weak, low-order 
perturbation theory works rather well. 



4 It ought be emphasized that all results with either the FD or the LS method represent converged eigenvalues, which means 
that contributions from folded diagrams to high order are included. 



G 



FIG. 1. Exact solutions (solid line) and the results obtained with the FD and LS methods with various approximations 
to the Q-box. FD-2nd and LS-2nd (dashed line) indicate that the FD and LS methods were used with a second-order Q-box, 
while FD-4th and LS-4th are the results obtained with a fourth-order Q-box. Similarly, LS-lOth and LS-50th are the LS results 
with a tenth- and fiftieth-order Q-box, respectively. A starting energy 0.5 was chosen in all calculations. 

For the LS method we note that with a low-order Q-box we are not able to reproduce the lowest eigenstate, even 
with an adequately chosen starting energy, as stated by the convergence criterion of the LS method. We actually 
have to go as far as to a Q-box of fiftieth order before we get a result close to the lowest eigenvalue. However, if the 
order of the Q-box is further increased (« 150 with the above parameters), even the LS method diverges. Thus, in 
summary, the general convergence properties of the LS method discussed in the introduction are demonstrated only 
for an exact Q-box. If the Q-box is approximated to a given order in the interaction, the convergence criterion cannot 
be applied. 

The implications of these results for the more realistic cases will be discussed in the next section. 

III. REALISTIC EFFECTIVE INTERACTIONS 

In calculations of effective interactions for 18 O, one normally chooses a model space consisting of a closed 16 O core 
with two nucleons in the single-particle states of the sci-shell. This model space consists of two-particle states only, 
and for J = we have three + states. However, from model calculations of energies and electromagnetic properties 
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and experimental data, the belief is that the 0^" state is predominantly of a four-particle-two-hole nature, while the 
ground state and the state are classified as two-particle states || . 

Before we proceed to study the convergence behavior of the folded-diagram and Lee-Suzuki expansions, let us first 
show how one can obtain useful information about the structure of the eigenstates from the slopes of the eigenvalues 
as functions of the starting energy lo. Using the property P + Q = 1, we can rewrite the energy dependent eigenvalue 
problem as 

/ PHP PHQ\( P*\ / P*\ . . 

^ QHP QHQ ) \ J a \ ) ' {lJ) 

where E is the exact eigenvalue. From eq. ( |l9[ ) we obtain 

E (Ptf |Ptf ) = (Ptf | PHP + PHQ 1— —QHP IP*} , (20) 

lo — QHQ 



with lo = E and 



such that 



mm) = (Py\HQ- ±—-pQH\P9), (21) 

(u - QHQ) 



( p *\ p *) = — ■ (22) 

This equation states that the slopes at the intersection between the r.h.s. of eq. (^0|) and the line E = lo provide 
information about the wave functions. The slopes 4^ should always satisfy 

au) 

in order to have (P\f r |P\I r ) < 1. Using H = Hq + Hi we can rewrite the r.h.s. of eq. ( p0| ) as 

PHP |*) = PH P |*) + PQP |*) = EP |*) , (23) 
where we have introduced the Q-box defined as 

PQP = PH 1 P + PH 1 Q ]——QH 1 P. (24) 

LO — Qrl Q 

How to evaluate the Q-box is however an open question. Here we will expand the denominator of eq. ( pi| ) such that 
we obtain a perturbative expansion for the Q _ box which reads 

PQP = PH X P + P f Hi — Q—Hi + Hi— 0— Hi— 0— Hi + . . ) P. (25) 
\ lj — Hq ui — Ho lo — Ho J 

The Q-box is made up of non-folded diagrams which are irreducible and valence linked. A diagram is said to be 
irreducible if between each pair of vertices there is at least one hole state or a particle state outside the model space. 
In a valence-linked diagram the interactions are linked (via fermion lines) to at least one valence line. Note that a 
valence-linked diagram can be either connected (consisting of a single piece) or disconnected, which in the case of a 
two-body interaction means that we have two one-body diagrams. In the FD and LS expansions, where we include 
folded diagrams as well, the disconnected diagrams are found to cancel out || . 

We evaluate the Q-box by including all two-body diagrams through third order in the reaction matrix G, as described 
m ref. §. The one-body diagrams plus the unperturbed term Hq are approximated by the experimental single-particle 
energies. Further, since we do not evaluate diagrams with Hartree-Fock insertions, we have no disconnected diagrams 
through third order in the interaction in our definition of the Q-box. The G-matrix replaces the interaction term Hi 
through Hi = G — U , U being an adequately chosen auxiliary potential, and is obtained through the solution of the 
Bethe-Goldstone equation 

G(lo) = V + V Q tt G{lo), (26) 
lo - Ho 
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where V is the free nucleon-nucleon interaction. Note that the Pauli operator Q may differ from the definition of Q 
in the above perturbative expansions. Eq. (^6|) is solved by using the so-called double-partition scheme described in 
ref. |[o]]. For V, we take the parameters of the Bonn B potential as defined in table A. 2 of ref. |ll[] with the Pauli 
operator Q defined so as to prevent scattering into states with one particle in the Os or Op shells or two particles in 
the IsOd or lpOf shells. A harmonic oscillator single-particle basis was employed with an oscillator parameter of 1.72 
fm. 



FIG. 2. The eigenvalues of eq. (20) for the three 0+ states of 18 0. 

In fig. 2 we show the graphical solutions for the 0^~, 0^" and 0^" states. As can be seen from fig. 2, the two lowest 
eigenstates exhibit negative derivates in the region near the graphical solutions, —12.2 and —8.3 MeV for the Of and 
the 0^ states, respectively. Correspondingly, the model space overlaps of eq. ( j22| ) in the region —20 to —4 MeV for the 
starting energy is of the order w 0.93 — 0.95 for the 0^~ state and ~ 0.98 — 0.99 for the 0^" state. We do not attempt 
to calculate the corresponding numbers near the crossing for the 0g state, since this would require the eigenvalue to 
be calculated at positive starting energies, which is not possible since the G-matrix is calculated at negative energies 
only. We are therefore not able to make any statement about the wave function overlaps for this state. In our further 
discussion we will therefore omit this state. 

For the 0^ state, the derivative becomes positive first at —28 MeV (4^ = 0.285), and the eigenvalue diverges as we 
approach the poles u> — QHqQ. However, it is interesting to see that within our approximation to the Q-box, the P- 
space overlap is close to one for the starting energies of interest, i.e. near the graphical solutions. This means that the 
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mixing from Q space states is rather weak within our perturbative scheme. As discussed above, the experimental 
state is expected to be of four-particle- two-hole structure (4p2h). Although we include the second-order 4p2h diagram 
in our definition of the Q-box, this contribution is not sufficient to represent the structure of the 4p2h correlations of 
the O2 state. The 4p2h diagram is also the only diagram in our third-order Q-box which may represent intruder state 
configurations. The omission of the latter diagram in our calculations does not change the qualitative aspect of the 
above wave function overlaps, indicating that we have a rather small Q-space component in conventional calculations 
of the effective interaction in the region of the graphical solutions. We may then summarize this to say that for a 
perturbative expansion up to third order in the G-matrix, the Q-space component is weak and one should identify 
the second eigenvalue of the model calculation with the 0^" experimental state. 

From the analysis in the previous section and the above results, we expect that the differences between the FD 
and the LS method, i.e. when comparing the converged results, should be small in the energy domain of interest. In 
particular, the P-space overlaps for the 0+ and the 0^ should be close to unity in both schemes. In fig. 3 we plot the 
eigenvalues as functions of the starting energy for the FD and LS methodsFl. 



FIG. 3. The eigenvalues derived with the FD and LS methods for the states of O. 



5 The one-body contributions plus Ho are also here approximated with the experimental single-particle energies 
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We observe that for starting energies within the range of the experimental values, the FD and LS methods show 
only small differences. The results here should therefore be interpreted in the same way as those exhibited in fig. 1. 
In that figure we showed that for low-order Q-boxes, the FD and LS methods gave essentially the same results. This 
leads us in turn to state that the convergence criterion for the LS method applies only if we are able to define an 
approximately exact Q-box. In realistic calculations this is however not the case, although both series converge rather 
rapidly, almost irrespectively of the choice of Q-box. Thus, in calculations of the spectra for, e.g. 18 0, we cannot 
choose a starting energy close to the 0^" state, in order to reproduce that state with the LS method, since the Q-box 
does not incorporate important Q-space configurations. Rather, the LS and FD methods give similar results, and the 
second + state obtained from the calculation should be identified with the experimental 0^ state, since this is believed 
to be predominantly a two-particle state. 

Finally, an important subtlety arises. Let us assume that we are able to define an exact Q-box for the two-body 
interaction. With an appropriately chosen starting energy, we should then be able to reproduce the intruder state 
in e.g. 18 employing the LS method. This interaction is then conventionally used in shell-model calculations to 
obtain observables of systems with more than two valence nucleons. However, for such heavier systems, the degrees 
of freedom of the intruder state in e.g. 18 may not be present. Thus, this two-body interaction obtained with the 
LS method should not be used in calculations of spectra of other sd-shell nuclei. With the LS method, we should 
therefore, in principle, recalculate the n-body interaction for the actual nucleus with n valence particles. This should 
be contrasted to the effective interaction obtained with the FD method. The latter reproduces only those states which 
have the largest P-space overlaps. If we restrict the attention to the degrees of freedom represented by the model 
space, we can use this two-body interaction in heavier systems as well. 



IV. CONCLUSIONS 



We have critically reviewed the convergence properties of the folded-diagram (FD) expansion of Kuo and co-workers 
[^| and the iterative scheme of Lee and Suzuki The behavior of these perturbative expansions was first studied 
within the framework of an exactly solvable model. Next, we investigated the convergence properties by applying the 
above two methods to a test nucleus, 18 0. Our conclusions are: 

• The exactly solvable model allows us to define an exact Q-box. With an exact Q-box we were able to demonstrate 
that, if we choose a starting energy close to a given eigenvalue, either P- or Q-space eigenvalues, the LS expansion 
converges to that eigenvalue. The FD method converges always to that eigenvalue which has the largest P-space 
overlap. 

• Still within the exactly solvable model, we have shown that if one chooses a Q-box of low order in the interaction, 
the FD and LS methods give almost similar results, irrespectively of the choice of starting energy. At low orders 
in the interaction we are not able to reproduce the Q-space state with the LS method, even with an appropriately 
chosen starting energy. The FD method (if it converges) always gives a result close to the eigenvalue with the 
largest P-space overlap. With increasing order of the interaction, the Q-box incorporates more Q-space degrees 
of freedom and the LS and FD methods give unstable results. 

• In realistic calculations, like those for 18 shown here, we can only define an approximate Q-box, to third order 
in the interaction in the present case. It is then seen that the FD and LS methods yield almost the same results, 
in agreement with the findings obtained from the simple model. Since the FD method converges to those states 
which have the largest P-space overlap, the structure of the LS wave function should therefore also have a 
large P-space overlap. With an approximative Q-box, we can therefore not choose a starting energy close to an 
intruder state in order to reproduce this state with the LS method. Thus, of the two methods, the FD method is 
the preferable one in nuclear structure calculations since it converges to those eigenstates which have the largest 
P-space overlap, a desirable property since we in general can only employ a rather small shell-model space. 



[1] M. Hjorth- Jensen, E. Osnes and H. Miither, Ann. of Phys. 213 (1992) 102 

[2] M. Hjorth- Jensen, T. Engeland, A. Holt and E. Osnes, Phys. Reports, in press 

[3] D.C. Zheng, J. P. Vary and B.R. Barrett, Nucl. Phys. A560 (1993) 211 



11 



[4] K. Suzuki, R. Okamoto, P.J. Ellis, J. Hao, Z. Li and T.T.S. Kuo, Phys. Lett. B308 (1993) 1 

[5] T.T.S. Kuo and E. Osnes, Folded-Diagram Theory of the Effective Interaction in Atomic Nuclei, Springer Lecture Notes 

in Physics, (Springer, Berlin, 1990) Vol. 364 
[6] K. Suzuki and S.Y. Lee, Prog. Theor. Phys. 64, 2091 (1980) 

[7] T.H. Schucan and H.A. Weidenmiiller, Ann. of Phys. 73 (1972) 108; 76 (1973) 483 

[8] P.J. Ellis and E. Osnes, Rev. Mod. Phys. 49 (1977) 777; P.J. Ellis and E. Osnes, Phys. Lett. B45 (1973) 425 

[9] P.J. Ellis and T. Engeland, Nucl. Phys. A144 (1970) 161 
[10] E.M. Krenciglowa, C.L. Kung, T.T.S. Kuo and E. Osnes, Ann. of Phys. 101, 154 (1976) 
[11] R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989) 



12 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9310013vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9310013vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9310013vl 



